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Abstract. Some models for quantum gravity (QG) violate Lorentz invariance and predict an en- 
ergy dependence of the speed of light, leading to a dispersion of high-energy gamma-ray signals 
that travel over cosmological distances. Limits on the dispersion from short-duration substructures 
observed in gamma-rays emitted by gamma-ray bursts (GRBs) at cosmological distances have pro- 
vided interesting bounds on Lorentz invariance violation (LIV). Recent observations of unprece- 
dentedly fast flares in the very-high energy gamma-ray emission of the active galactic nuclei (AGNs) 
Mkn 501 in 2005 and PKS 2155-304 in 2006 resulted in the most constraining Hmits on LIV from 
light-travel observations, approaching the Planck mass scale, at which QG effects are assumed to 
become important. I review the current status of LIV searches using GRBs and AGN flare events, 
and discuss limitations of light-travel time analyses and prospects for future instruments in the 
gamma-ray domain. 
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INTRODUCTION 

Currently we have two extremely successful fundamental theories at hand in physics: 
general relativity, describing space, time and gravitation; and quantum theory, describing 
electroweak and strong interactions by means of quantum field theories. The interface of 
both theories was not explored for a long time, because there have been no experiments 
that could probe regimes in which both quantum and gravitational effects become 
important. Still, not only from a philosophical point of view, physicists aspire a unified 
theory of both, which is usually called quantum gravity (QG). In standard relativistic 
quantum field theory (QFT), space-time is considered a fixed arena in which physical 
processes take place. Quite in contrast, QG theories require drastic modifications of 
space-time when structures approach the QG mass scale [1, 2, 3], often identified with 
the Planck mass Mpianck = \/^cjG^ = 1.22x 10^^ GeV. However, the energies at which 
QG effects lead to deviations from conventional quantum mechanics might lie well 
below the QG energy scale. As a non-renormalizable interaction gravity may leave 
distinctive imprints also at energies much below the preferred QG scale, if violating 
any fundamental symmetry. Thus, we might spot QG effects not only on extreme energy 
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scales, but use GeV/TeV photons as "low-energy" probes. 

An energy dependence of the speed of light in vacuum may arise from photon 
propagation through a gravitational medium containing quantum fluctuations on dis- 
tance scales on the order of the Planck length /pianck = s/fiGj^/c^ ~ 10~^^cm and on 
timescales ^ ^pianck- While being smooth at large distances, space-time at short dis- 
tances of the order of /pianck might show a complex, foamy, structure due to quantum 
fluctuations [4, 5, 6]. Such Planck-size topological defects, viz. black holes with micro- 
scopic event horizons, appearing and evaporating spontaneously, might lead to quantum 
decoherence and Lorentz invariance (LI) violation. There are numerous formal mathe- 
matical approaches to QG. Among the most popular ones are loop QG [7, 3], which is 
the canonical non-perturbative quantization of general relativity, and superstring theory, 
which originally has been developed to describe the behavior of hadrons. (See [8] for 
a review on QG theories). In Liouville string theory [9], modified dispersion relations 
(MDR) are enabled by higher-order derivatives in effective Maxwell and Dirac equa- 
tions. 

It has been pointed out by Amelino-Camelia et al. [1] that different approaches to QG 
lead to a similar description of first-order effects of such a time dispersion: 

Ar^^-^- (1) 

£^QG C 

where A? is the time delay relative to the standard speed of light, ^ is a model-dependent 
factor of the order 1, £■ is the energy of the observed radiation, and L is the distance 
traversed by the radiation. Thus, LI violation (LIV) is expected as a generic signature 
of approaches to QG. A rather general review on the status of LIV tests is given by 
Mattingly [10]. 

A number of predicted new phenomena lead to testable astrophysical effects. Dis- 
persion [1] and vacuum birefringence [3] (leading to polarization-dependent LIV) are 

purely kinematical effects and require only an MDR and a standard definition of the 
group velocity. Below, we will particularly elaborate on astrophysical time-of-flight tests 
of the dispersion relation of photons. Other effects, like anomalous threshold reactions, 
threshold shifts in standard reactions, or reactions affected by "speed limits" (e.g. syn- 
chrotron radiation), need additional assumptions on energy/momentum conservation and 
dynamics, or on effective QFT, and will not be discussed here. Note that a non-constant 
speed of light may also occur with maintained Lorentz invariance, e.g., in doubly-special 
relativity [11]; this underlines the necessity for astrophysical tests of LIV (LIV tests on 
Earth would yield negative results!). 

In the remaining sections, we will review astrophysical probes of QG with high energy 
y-rays. 



PROBING QUANTUM-GRAVITY WITH PHOTON 
TIME-OF-FLIGHT MEASUREMENTS 

When looking for QG phenomena, we expect deviations from QFT, presumably sup- 
pressed by some power of the Planck mass (or the QG mass scale). From a purely 
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phenomenological point of view [1], such effects can be treated using a perturbative 
expansion, assuming the involved energies E <C Mpiancki 




(2) 



Note that an explicit breaking of LI is expected at the Planck mass scale. A linearly- 
deformed dispersion relation arises, e.g., in noncritical Liouville string theory, while in 
critical string theory, only quadratic deviations are expected [12]. An MDR implies an 
energy-dependent speed of light. 



or, in other words, the vacuum acquires non-trivial optical properties, namely a refractive 
index v{E) = c/n{E). This translates into a delay, including cosmological effects, 

^ = «i^^«i^"o7W (4) 

To measure any such time delay, very fast transient astrophysical phenomena are re- 
quired, which provide a time stamp for the simultaneous emission of photons of dif- 
ferent energies (Fig. 1). The corresponding figure of merit for QG tests is given by the 
sensitivity to a given QG mass scale, 

Mqo = «^^ (5) 

The photon energy E acts as lever arm, implying that both an instrument for measuring 7- 
ray energies as high as possible, but also astrophysical sources that deliver these energies 
are advantageous. But also the time resolution of the instrument and the fineness of the 
structures in the y-ray signal, described by At, are decisive, and the source has to provide 
enough photon statistics per time as to allow for a binning as small as possible. L denotes 
the luminosity distance of the source. 



QUANTIFICATION OF TIME DELAYS 

In the X-ray to VHE y-ray domain, photon data are characterized by arrival directions, 
energies, and arrival times of X-rays, gamma rays, or gamma-ray candidates. A very 
basic method of assessing any time delay of very high energy photons, motivated also 
by the sparse statistics at these energies, relies on finding peaks in the cross-correlation 
function of light curves of different photon energy bandpasses [13, 14, 15]. Note, 
however, that such analyses do not take into account the full available information 
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FIGURE 1. In photon time-of-flight measurements, it is assumed that the energy dispersion at the 
source is known or neghgible. Any additional dispersion can then be attributed to the dispersive prop- 
erties of the vacuum, and this dispersion can be measured by the observer. Original drawing courtesy 
A. MorselH. 



each individual photon carries [16]. Therefore, more appropriate methods have been 
developed, as summarized in the following. 

• In the energy-cost-function method [16], it is assumed that traveling through any 
dispersive medium will increase the spread of the photon arrival time distribution.^ 
Thus, by minimizing the spread of an observed, dispersive-broadened pulse by an 
"energy cost function", the original pulse and the dispersion can be reconstructed 
[16, 17]. 

• A log-likelihood fit, describing the flare peak form Oe (which might be assumed, 
e.g., as Gaussian pulse), the instrumental energy smearing G, the energy spectrum 
T{Es) of the flare, and the observed photon arrival times Fs{ts) can be applied to the 
unbinned data [18, 19], e.g., 

diP n 

-— = / r(£,)G(£ aE{E,,))F,{t,)dE,. (6) 
at at Jo 

• Cross-correlation of oversampled light curves [20]. This "modified cross- 
correlation function" (MCCF) method allows for resolving features (as time 
delays) in the light curves well below the duration of the flux bins. 

• A wavelet analysis suppresses noise in the data and identifies "genuine variation 
points" [21], which can be attributed to flares and can be used as time stamps in 
different energy bandpasses. This method is particularly useful if the observational 
data are noisy and the (multiple or complex) flare structure is not very significant. 



Note that the dispersion needs not be source-extrinsic, thus even dispersive regions inside a y-ray source 
can be modeled using the ECF approach. 
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PRESENT REACH USING HIGH-ENERGY GAMMA-RAYS 

Pulsars have very distinct periodic light-curve forms, and due to the periodicity, these 
light curves can be sampled very accurately. In observations of photons up to E 2GeV 
from the Crab pulsar, at a distance of 2.0 ±0.5 kpc [22], using the sub-millisecond 
timing of EGRET, A? = 10^ s was achieved. With these numbers, a sensitivity to Mqq « 
lO^^GeV can be achieved. For active galactic nuclei (AGN; typical distances: in excess 
of 100 Mpc = lO^'^cm), Imaging Air Cerenkov Telescopes (lACT [23, 24, 25]) currently 
may observe bright y-ray flares with energies up to around Ej ^ lOTeV with sufficient 
statistics and a time resolution (given by the AGN flare timescales) of ~ 50 — 100s. This 
enables a sensitivity up to the Planck energy scale, Mqq ^ lO^^GeV. In the absence of 
a linear term in the MDR, the quadratic reach would be around Mqq = 6 x lO^GeV. 
Gamma-ray bursts, at distances beyond 7000 Mpc, will be routinely observed by Fermi 
[26] at energies up to i^(GeV), which will yield a sensitivity of Mqq 10^^ GeV. Note, 
however, that current GRB studies are restricted to Mqq — 10^^ GeV, as until recently 
only AE lOOkeV was available {Mqq = 10^ GeV in the quadratic case). 

ASTROPHYSICAL VHE GAMMA-RAY SOURCES 
Pulsed emission from the Crab nebula 

The Crab pulsar has a pulsation period of 33. 18 ms. As its pulsations are well aligned 
in time from radio through X-ray wavelengths, it seems likely that photons of different 
energies are produced nearly simultaneously. Using the sub-ms timing of the EGRET 
experiment, it could be shown [27] that > 2 GeV photons trail those at 70—100 MeV 
by no more than 0.35 ms (95% c.L), implying Mqg > 1.8 x 10^^ GeV. Recently, pulsed 
emission from the Crab was also detected at energies up to 60 GeV [28, 29]: About 
8,500 7-ray events form a signal at the 6.4-a significance level. This observation might 
yield a lower limit for quantum gravity effects of Mqg 1 .2 x 10^^ GeV [30], exceeding 
the limit given by the EGRET-Crab analysis. 

Gamma-Ray Bursts 

Gamma-ray bursts (GRBs; [31]) have rather complex individual time profiles. Once a 
relativistic fireball is created in the gigantic explosions thought to be the origin of GRBs, 
the physics ought to be insensitive to the nature of the progenitor object. GRBs have 
been found at redshifts exceeding z = 5 (7% of the presently known GRBs; the average 
redshift is z = 2.8). Thus GRBs seem to be a natural proxy to test QG predictions. The 
key issue is to distinguish possible QG signatures from intrinsic delays, particularly as 
delays in GRBs have been observed. QG effects should increase predictably with the 
distance of the individual GRB. The accessible energies, however, have been smaller 
than or at most some GeV; lACT, which are sensitive in the above 30 GeV region, 
are still awaiting the first GRB detection [32, 33]. An additional complication arises 
from the fact that establishing the distance of GRBs is not trivial and works only for a 
fraction of the observed GRBs. The strongest QG constraint from an individual GRB is 
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that from GRB 930229 [34], where during a rise time of 220 ± 30/is the dispersion of 
photons was measured to be Ar < 25 ms at energies Afi = 30keV — 80MeV, yielding 
Mqg > 8.3 X IQi^GeV. The caveat of this observation is that GRB 930229 has no 
measured redshift(!), thus with an assumed peak flux of 10^^ photons s~^ and the peak 
width-redshift relation [35] a distance D = 260 Mpc was inferred. For GRB 051221A, 
which occurred in a galaxy at z = 0.55, no time delay of photons between 15 and 350 
keV was observed [36], leading to a limit of Mqq > 6.6 x lO^^GeV. 

A more robust method is to study large samples of GRB. In [37], a total of 36 GRBs 
observed by BATSE (9 bursts), HETE-2 (15 bursts), and Swift (11 bursts) up to z = 6.3 
was used. The energy difference of photons in that sample is around 200 keV, and 
to account for possible source effects, a "universal source uncertainty" of 54 ms was 
accounted for. A limit of Mqq > 0.9 x lO^^GeV was deduced. 

Fermi Large-Area Tracker (LAT) and Gamma-Ray Bursts. In LAT, around 250 
bursts per year are expected [38] to be detected, and about 50 of them will survive a 
high-energy cut, i.e., contain > 550 MeV photons. LAT itself should detect around 15 
bursts without Gamma-ray Burst Monitor triggers. Kuehn et al. [39] simulated the LAT 
response to GRB photons, focusing on the prompt GRB emission. Due to the wide field 
of view, LAT in contrast to the TeV instruments is quite unique in detecting prompt 
GRB emission. For the simulations, spectral indices have been taken from BATSE and 
extrapolated to LAT energies. Even with a binned analysis, a simultaneous fit to multiple 
GRBs (at z = 1) was shown to have a sensitivity at 10^^ GeV. An alternative method [17] 
also has this sensitivity. Recently, the LAT collaboration reported the first detection of 
high-energy photons from three GRBs [40], and used the strongest burst, GRB 080916C, 
to derive a lower QG energy scale limit of Mqq > (1.5 ±0.2) x 10^^ GeV. 

Flares from Active Galactic Nuclei 

The study of high energy (E > 100 MeV) y-ray emission from active galactic nuclei 
(AGN) is one of the major goals of space-borne and ground-based y-ray astronomy. 

The vast majority of the currently known VHE y-ray emitting AGNs [41] are blazars 
[42],^ a subclass characterized by high variability over the whole electromagnetic spec- 
trum, a characteristic peak in their spectral energy distributions (SEDs) from nonthermal 
synchrotron emission, and a second peak at y-ray energies. In synchrotron- self -Compton 
(SSC) models it is assumed that the observed y-ray peak is due to the inverse-Compton 
(IC) emission from the accelerated electrons up-scattering previously produced syn- 
chrotron photons to high energies [43]. In hadronic models, instead, interactions of 
a highly relativistic jet outflow with ambient matter, proton-induced cascades, or syn- 
chrotron radiation off protons, are responsible for the high-energy photons. 

The duty cycle of Fermi enables a quasi-continuous monitoring of blazar flaring ac- 
tivity (also monitoring programs by lACTs are in place, [44]), while the high sensitivity 
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of current lACT is advantageous for understanding the particle acceleration mechanisms 
in flares, the origin of the high-energy y-rays, as well as the relations between photons 
of different energies (from radio to VHE; particularly between the two emission popu- 
lations [45]). 

The Whipple QG Limit from a Mkn 421 Flare 

The first rapid flare from a TeV blazar was observed by the Whipple telescope from 
Markarian (Mkn) 421 (z = 0.030) on 1996 May 15 [46]. The peak of the flare is 
concentrated in one 280 second bin, and as there were no > 2 TeV /-rays observed 
outside this bin, it was concluded that the y-rays < 1 TeV and > 2 TeV are in step with 
>95% probability [14]. This observation translates in a lower limit Mqg > 4 x 10^^ GeV. 

The July-2005 Flares of Mkn 501 

MAGIC [23] recorded fluxes exceeding four times the Crab-nebula flux from Mkn 
501 (z = 0.034) in 2005, and revealed rapid flux changes with doubling times as short 
as 3 minutes or less [15]. For the first time, short (~ 20 minute) VHE y-ray flares with 
a resolved time structure could be used for detailed studies of particle acceleration and 
cooling timescales. The two flares behaved differently (Fig. 2): While the 2005 June 
30 flare is only visible in 250GeV-1.2TeV, the 2005 July 9 flare is apparent in all 
energy bands (120 GeV to beyond 1.2 TeV). Additionally, at a zero-delay probability of 
P = 0.026, a marginal time delay At = ZiE with T; = (0.030 ±0.012)sGeV^ towards 
higher energies was found using two independent analyses, both exploiting the full 
statistical power of the dataset (see [16, 19] for details). For a quadratic effect (vanishing 
first- order term in eq. 2), At = XqE'^ with Xq = (3.71 ±2.57) x lO^^sGeV"^ obtained 
correspondingly. Several explanations for this delay have been considered up to now: 

1. Electrons inside the emission region moving with constant Doppler factor are 
gradually accelerated to energies that enable them to produce corresponding 7 rays 
[15]. 

2. The 7-ray emission has been captured in the initial acceleration phase of the rel- 
ativistic blob in the jet, which at any point in time radiates up to highest 7-ray 
energies possible [47]. 

3. An one-zone SSC model, which invokes a brief episode of increased particle 
injection at low energies [48]. 

4. When assuming a simultaneous emission of the 7-rays (of different energies) at 
the source at z = 0.034, an observed time delay can be converted in QG mass 
scales Mqgi = 1.445 x lO^^s/l/ and Mqg2 = 1-222 x 10^(s/t^)1/2, respectively. 
As the measured time delay is marginal, they translate in a lower limit of Mqgi > 
0.21 x 10^^ GeV (95% c.l.) for the linear case and Mqg2 > 0.26 x lO^^ GeV (95% 
c.l.) for a quadratic dependence on energy [16]. Results for the 2005 June 30 exhibit 
a similar sensitivity, but this flare is not very significant. Both limits increase further 
if any delay towards higher energies in the source itself is present. 
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FIGURE 2. Light curves for the observations on 2005 June 30 (left panels) and July 9 (right panels) 
[15] with a time binning of 4 minutes, and separated in different energy bands ranging from 150 GeV to 
>1.2 TeV. For comparison, the Crab nebula emission level is shown. The data are divided in into 'stable' 
and 'variable' emission. The probabilities given in the left panels are for fits with constant lines in the 
'variable' region. Only the constant-line fit in the energy range 0.25-0.6 TeV is not satisfactory. In the 
right panels, the 'variable' regions were fitted with a flare model. A more quantitative photon-by-photon 
analysis of these data was performed in [16]. 



The 2006 Giant PKS 2155-304 Flare 

In the night of 2006 July 28, the H.E.S.S. array [24] observed a giant (> 10 Crab) 
outburst of the blazar PKS 2155-304, almost four times as distant as IVIkn 421 and IVIkn 
501 (z = 0. 116). Also in this case, fast variability on the order of 3 minutes and an energy 
coverage up to a few TeV is given [49]. For searching for time delays, a IVICCF method 
[20] and a wavelet method were employed, both yielding a zero time delay within errors, 
and thus T, > 0.073 sGeV^^ and = 45 x lO^^sGeV"^, respectively (95% c.L). This 
translates into lower limits of Mqgi > 0.52 x 10^^ GeV and Mqqi > 0.014 x 10'° GeV 



Source-Intrinsic Lags and the Blazar Ensemble 2008 

IVlost acceleration and cooling mechanisms scale in energy, although there is no 
definite need for a strong energy dependence. In X-ray observations, however, soft and 
hard lags have been observed (most recently, e.g., by [50]), and lead to characteristic 
loops in flux-spectral hardness plots [51]. Those diagnostic plots may also be used 
in the VHE regime [52]. Intrinsic lags may strongly depend on the individual blazar 
and even on the individual flare, as the three AGN outburts used for QG searches have 
proven. With our current understanding of blazar physics, intrinsic delays cannot be 
corrected for individually. They will, however, not scale with the source distance, so 
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FIGURE 3. The PKS 2155-304 flare [49]. Black points show flie Ught curves measured between 
280 — 800 GeV (upper panel) and >800 GeV (lower panel), binned in two-minute time intervals. Gray 
points show the oversampled light curve using the MCCF method [20], for which the two-minute bins are 
shifted in units of five seconds. The inlay in the upper panel illustrates this in a zoom, where the horizontal 
error bar shows the duration of the bin in the original light curve. 
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FIGURE 4. The blazar flare ensemble 2008: Three limits on the energy dispersion as measured in flares 
of Mkn 421, Mkn 501, and PKS 2155-304. If die marginal time delay in Mkn 501 is extrapolated, it is 
found to be incompatible with the null result for PKS 2155-304. 



robust conclusions on QG effects ultimately require the study of samples of very fast 
flaring objects at different redshifts. These are expected to be observed by the current and 
next- generation [53] lACTs and Fermi thanks to the high sensitivity of these instruments. 
Note that multi-wavelength observations of flares are complementary to QG searches 
and essential, as they will allow to study possible source-intrinsic mechanism and effects 
and enhance the understanding of AGN source physics. Any such understanding will 
enhance QG limits. 

Note further that obviously source-instrinsic and QG effects may cancel out, leading 
to an observation of a null effect. This scenario has been dubbed "conspiracy of Nature" 
[14] and while it is still conceivable that both effects do exhibit the same timescale, but 
the opposite sign, also here a sample of blazar flares from blazars at different redshifts 
will help to separate an eventual QG effect from blazar physics. Fig. 4 shows the current 
blazar flare ensemble observed in VHE y-rays. 
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THE FUTURE: PROSPECTS AND LIMITATIONS 



High energy observations of fast pulses enable measurements of the energy dependence 
of the speed of light. Such a possible LIV is predicted by many QG models. An energy- 
dependent lag of TeV 7-rays on the 2.5-a level has been observed in a flare of the blazar 
Mkn 501 . The constraints derived from the ~ 4 times more distant blazar PKS 2 1 55-304, 
however, are consistent with a null result, as are robust ensemble analyses from GRBs 
at a lower Mqq sensitivity. All in all, VHE observations of blazars constitute the most 
stringent limits on a linear (PKS 2155-304) and quadratic (Mkn 501) energy dependence 
in the speed of light. Future improvements in VHE searches may be expected from: 

At - Even the high time-resolution lACTs have not yet explored the shortest blazar 
variability timescales, although high photon statistics are available; e.g., more than 
10,000 photons in the PKS 2155-304 flare. Thus, only "brighter flares" or even 
more sensitive instruments will help to improve on At. Observing faster variations, 
or flares/subflares with shorter timescales (say 30 seconds instead of currently 3 
minutes), might increase the QG sensitivity by about x5. GRBs with their fine 
time structures are interesting to observe, particularly with Fermi. 

AE - Gaining access to higher energy y-rays enlarges the lever arm in the QG sen- 
sitivity. However, the bright VHE blazars usually exhibit rather steep spectra, 
~ E~'^-^...E~^-^ [41], thus it is rather unlikely to gain an order of magnitude more 
sensitivity, perhaps a factor x3 is feasible. 
L- It it difficult to improve much using the known bright, close-by blazars, as we have 
not yet seen flares exceeding the 15-Crab level, even not during the giant 1997 flare 
of Mkn 501 [54]; More distant blazars, however, seem to exhibit steep spectra (e.g. 
3C 279 [55, 23]), which are expected due to y-ray attenuation by the extragalactic 
background light [56, 57, 55, 58]. A major flare in, e.g., 3C 279, extending the 
distance from z = 0.12 to 0.54, corresponds to a gain by x5 in sensitivity. As 
discussed above, for robust ensemble studies, we need to find flares at various 
redshifts; thus, improved lACT facilities [53] to reach lower energies than the 
currently reachable (fa 30 GeV by MAGIC; « 100 GeV by VERITAS/H.E.S.S.). 

The overall increase to be hoped for with future lACT observations is thus about x70. 

Fermi will regularly detect distant (z > 1) GRBs, thus collecting statistics and be able 
to build a robust GRB sample. This will increase the sensitivity for such a robust GRB 
study by about x20. 
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